function [vhx] = VHx(xh,xcd,VE,VHy)
vhx=sign(xcd-xh)*power(power(VE,2)-power(VHy,2),1/2);
end

